Observation of Low-Frequency Raman Modes and Electronic Excitations In 

Atomically Thin MoS 2 Crystals 
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Atomically thin M0S2 crystals have been recognized as a quasi-2D semiconductor with remark- 
able physics properties. This letter reports our Raman scattering measurements on multilayer and 
monolayer M0S2, especially in the low-frequency range (<50 cm -1 ). We find two low-frequency 
Raman modes with contrasting thickness dependence. With increasing the number of M0S2 layers, 
one shows a significant increase in frequency while the other decreases following a 1/N (N denotes 
layer-number) trend. With the aid of first- principle calculations we assign the former as the shear 
mode E\ g and the latter as the compression vibrational mode. The opposite evolution of the two 
modes with thickness demonstrates novel vibrational modes in atomically thin crystal as well as 
a new and more precise way to characterize thickness of atomically thin M0S2 films. In addition, 
we observe a broad feature around 38 cm -1 ( 5 meV) which is visible only under near-resonance 
excitation and pinned at the fixed energy independent of thickness. We interpret the feature as an 
electronic Raman scattering associated with the spin-orbit coupling induced splitting in conduction 
band at K points in their Brillouin zone. 

PACS numbers: 63.22.-m, 78.30.-j, 71. 15. Mb, 71.35.-y 



In the hunt for ultimately thin electronic devices, 
atomically thin molybdenum disulfide (M0S2) layers have 
been recognized as a promising material after the success 
in graphene. Bulk M0S2 crystal has a structure of S-M- 
S covalently bonded hexagonal quasi-2D network weakly 
staked by weak Van der Waals forces as sketched in Fig- 
ure la. M0S2 monolayer, the elementary unit to form 
ultrathin films by weak stacking is stable under ambient 
conditions and is easily fabricated either by mechanical 
or liquid-phase exfoliation-^—. As M0S2 thins to a few 
atomic layers, there emerge many exotic physical prop- 
erties absent in bulk crystal. For example M0S2 thin film 
undergoes a transition from an indirect-gap semiconduc- 
tor with a gap of 1.2 eV in bulk form to a direct- gap 
one with a band gap of 1.9 eV in monolayers -^. These 
distinct physical properties make atomically thin M0S2 
crystal attractive for diverse applications such as elec- 
tronic devices, spintronics, photovoltaic, energy storage, 
etc. 

Raman scattering has been proven to be valuable as a 
nondestructive and versatile method of structural char- 
acterization and as a method for studying electronic and 
vibrational properties in quasi-2D materials. There are 
four Raman-active vibrational modes in bulk M0S2 crys- 
tal namely Ei g (in-plane optical vibration of S atoms), 
Ai 9 (out-of-plane optical vibration of S atoms), E\ g (in- 
plane optical vibration of Mo-S bond) and E\ g modes 
(in-plane optical vibration of rigid atomic bond) 6 . The 
Ei g (in-plane optical vibration of S atoms) of M0S2 thin 
films is negligible in the back scattering geometry which is 
widely used in micro-optic spectroscopic study in nanos- 
tructures, owing to the forbidden selection rule from the 
symmetry point view. Ryu et at investigated the in- 
plane E\ g mode and out-of-plane h\ g mode around 400 
cm -1 in few-layer M0S2 and found that the two modes 



show opposite frequency dependence on the film thick- 
ness, E\ g softening and K\ g stiffening with thickness pre- 
sumably due to increased dielectric screening, stacking- 
induced changes in intralayer bonding and interlayer cou- 
pling. The frequency difference of these two characteris- 
tic E\ and Ai g modes ranging from 19 cm -1 in mono- 
layer to 26 cm -1 in bulk crystal are well accepted as a 
thickness monitor for atomically thin M0S2 flakes. The 
E\ g mode (in-plane optical vibration of rigidly bonded 
atoms) also known as shear mode originates from the 
rigid in-plane vibration against adjacent layers and seem- 
ingly exist in many layered structures. The energy of the 
E\ g mode is very low&, analogously 30-40 cm -1 in few- 
layers grapheneS-, and is hardly observable in many com- 
mercial apparatus due to the limited rejection against 
Rayleigh scattering. 

Here we report our observation of Raman scattering 
on few-layer M0S2 in the low frequency range. We find 
two low-frequency Raman modes (l~6meV) with con- 
trasting thickness dependence. With increasing thick- 
ness of M0S2 films, one shows a significant increase in 
frequency while the other decreases following a 1/N (N 
denotes the number of unit layers) trend. With the aid 
of first-principle calculations we assign the former as the 
shear mode E\ and the latter as the compression vibra- 
tional mode. The opposite evolution of the two modes 
with thickness demonstrates novel vibrational modes in 
atomically thin crystals as well as a new way to precisely 
determine the number of M0S2 atomic layers. In addi- 
tion, we observe a broad feature around 5 meV, which is 
visible only under near-resonance excitation and locates 
at the fixed position at various thicknesses. 

Ultrathin flakes were mechanically exfoliated from a 
natural M0S2 single crystal onto silicon wafers capped 
with a 300nm thick SiC>2 by a method analogous to that 
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FIG. 1: (Color online) (a) Optical images of ultrathin M0S2 flakes, (b) Schematics of Raman active modes, (c) Low-wavenumber 
Raman spectra of atomically thin M0S2 flakes with different thickness under 532nm excitation. The dashed line denotes the 
evolution of the weaker low-frequency mode, (d) The corresponding E 1 2 9 and Ai 9 modes. Their frequency difference is a 
measure of the number of layers. 



of producing graphene. Atomically thin M0S2 flakes were 
first visually screened with interference color through op- 
tical microscope. Typical optical images of M0S2 ultra- 
thin slabs are presented in Figure 1. The film thickness 
is confirmed by atomic force microscope and photo lumi- 
nescence spectra 4 * 5 -. The Raman scattering is collected 
with a confocal setup with a triple stage monochromator 
with excitation at 532nm and 633nm. 

Figure 1 presents the Raman spectra from samples 
with various thicknesses. The in-plane vibrational E\ g 
mode and the out-of-plane vibrational Ai g mode around 
400 enr 1 show a consistent trend as reported in Ref.7. 
The E\ g mode stiffens and the Ai g mode softens as the 
sample thins to multi- and monolayers. The frequency 
difference between the two mode shows a consistent indi- 
cation of the layer thickness as other means, particularly 
for monolayer, bilayer, trilayer and quadlayer (N<4). In 
the low-frequency range, two modes can be clearly seen 
at various thicknesses except monolayer at 532nm ex- 
citation. The mode with narrow bandwidth (^2cm _1 ) 
undergoes a blue shift with increasing thickness, from 22 
cm -1 in bilayer to 32 cm -1 in bulk form, while the mode 



with broad bandwidth undergoes an opposite evolution 
as indicated by the dash line in figure 1. 

To gain insight into the observed modes, we have stud- 
ied the zone center phonons of M0S2 from first-principles 
density functional theory (DFT) calculations as imple- 
mented in the Vienna ab-initio simulation package 1 ^—. 
The few-layer M0S2 were modeled as a lxl unit cell 
in 2~5-layer slabs, which were separated by a vacuum 
>13 A in the perpendicular direction. The projector 
augmented wave (PAW) metho d 14 ' 15 and the local den- 
sity approximation (LDA) for exchange-correlation po- 
tentials were adopted. The energy cutoff for the plane 
waves was 400 eV. The first Brillouin zone was sampled 
with a 12x12x4 k-point mesh for bulk and a 12x12x1 
mesh for few-layers. The Gaussian broadening with a 
width of 0.05 eV was used for Fermi surface smearing. 
All ions were fully optimized until forces were less than 
0.01 eV/A. The vibrational frequencies and displacement 
patterns of phonon modes were calculated at the equilib- 
rium structures using the dynamical matrix method as 
described in our previous publication 16 . It is worth not- 
ing that although the LDA method neglects the van der 
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FIG. 2: (Color online) Calculated phonon frequencies and 
displacement patterns for the shear modes of (a) bulk, (b) 
bilayer, (c) trilayer, (d) quadlayer, and (e) pentalayer M0S2 
at the Brillouin zone center obtained from DFT calculations. 
The Raman-active (R) and infrared-active (I) modes are la- 
beled. Atomic structures and displacement patterns prepared 
using XCRYSDENfi^ Modes for different layers are high- 
lighted with dashed boxes. 



Waals interaction between the layers, it produces reason- 
ably well the inter-layer phonons of MoS^. 

The calculated frequencies of E\ g (392.7 cm -1 ) and Ai ff 
(419.3 cm -1 ) modes at the Brillouin zone center in bulk 
form show good agreement with our experimental values 
(382.9 cm -1 and 408.2 cm -1 ), as well as with previous 
experiments- (382.0 cm -1 and 407.5 cm -1 ) and theo- 
retical results (387.8 cm -1 and 412.0 cm -1 )^. Here we 
concentrate on the in-plane shear modes below 50 cm -1 
which exhibit prominent thickness-dependence as shown 
in Figure 1. Figure 2 presents the calculated phonon 
frequencies and displacement patterns for the double- 
degenerate shear modes (E\ g mode) in bulk and few-layer 
M0S2. The calculated shear modes: 35.3 cm -1 in bi- 
layer(2L), 37.8 cm -1 in trilayer(3L), 38.8 cm -1 in quad- 
layer(4L) and 39.0 cm -1 pentalayer(5L) agree qualita- 
tively well with the corresponding observed sharp peaks 
around 30cm_i in Figure 1. The phonon frequencies of 
these in-plane shear modes increase from null in mono- 
layer monotonically with the layer thickness. They show 
similar thickness-dependence and displacement patterns 
as the shear modes in multilayer graphene 9 , illustrating 




N (Number of MoS 2 layers) 

FIG. 3: (Color online) (a) Thickness dependence of fre- 
quencies. The blue curve is the fitting with the formula 
u) = yg,^ ^|\A + cos ( j|0- ( D ) Thickness dependence of 
compression mode frequencies (see below) . The red and green 
curves are produced with the least-square fitting. 

their similar origin from the interlayer interactions. So we 
can safely assign the sharp peaks in Figure 1 to the shear 
modes (E^g mode as denoted in bulk) of adjacent lay- 
ers oscillation driven by Van de Waals interaction. Our 
assignments are consistent with previous reports in bulk 
and few-layer MoS^ 1 ^. 

As such a low-frequency interlayer shear mode is a ge- 
netic vibration for layered compounds, we examine the 
peak position of the shear modes with a linear-chain 

model 9 : u — ^ 7ic W^V^ + cos ( j^) where w, c, a, [i and 

N denote the peak position of the shear mode in cm , 
speed of light in cm/s, the interlayer force constant per 
unit area, unit layer mass per unit area and the number 
of layers, respectively. 

Figure 3a shows a perfect fitting with the interlayer 
force constant per unit area a = 2.9 x 10 19 N • m~ 3 and 
the single-layer mass per unit area \i = 30.3kg- A~ 2 . One 
can further estimate the shear modulus C44 — a xt with 
the spacing t between adjacent M0S2 layers is ~ 6.2 A. 
This immediately gives C44 ~ 17.9 GPa. The C44 value is 
in good agreement with the reported experimental shear 
modulus in M0S2 bulk crystal, from 15 to 19 GP a 20 ' 21 . 

Now we move to the other low-frequency Raman mode. 
To our knowledge, this broad mode has not been reported 
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FIG. 4: (Color online) Low-wavenumber Raman spectra of 
M0S2 thin Ems at different thicknessunder near-resonant ex- 
citation (633nm). The broad feature around 38 cm -1 is ob- 
served only in near-resonant excitation. 



before. As mentioned above, it cannot be assigned to 
any known Raman-active mode. It features a contrast- 
ing thickness-dependence against thei?f g mode, as shown 
in Figure 3b. Besides the thickness-dependence follows a 
perfect 1/N behavior rather than 1/yN. The signature 
of 1/N behavior has been observed in similar thin- film 
systems with a number of atomic layers grown by molec- 
ular beam epitaxy (MBE), such as Na/Cu(001)22. And 
it is well understood as open-ended standing wave mode 
vertical to atomic planes with a linear-chain force con- 
stant model, particularly including the coupling between 
atomic layers and substrate^. For layered compounds, in 
principle there should exist a compression mode besides 
the shear mode discussed above, which is a rigid-layer 
vibration vertical to layer. The longitudinal-like open 
standing wave will give an effective relative shift between 
adjacent layers if one looks at its vibration amplitude, as 
shown in Figure lb. That means it actually corresponds 
to the compression mode. The compression mode, com- 
bined with a longitudinal optical phonon mode, has been 
revealed in few-layer graphene through double-resonant 
intravalley scattering^. The mode strongly depends on 
the coupling between atomic layers and substrate. And 
the Raman intensity of the mode is much smaller than 
that of the shear mode. This may be the reason why it is 
not easy to be observed in few-layer graphene and M0S2. 

The direct observation of the shear mode and the com- 
pression mode provides a new way to more precisely char- 
acterize the thickness of ultrathin M0S2. A conventional 
method is to measure frequency difference between Ai ff 
andi?2g mode. Considering the contrasting evolution of 
the two rigid-layer low-frequency modes with thickness 
and the large relative shift up to 100%, the frequency 
difference between the shear mode and the compression 
mode would be a more precise measure of layer number, 
particularly for slabs with more than 3 unit layers. 



Another new finding in our Raman measurements is 
the broad feature centred at 38 cm -1 (~5meV) mea- 
sured with a cross-polarization configuration under a 
near-resonant excitation, as shown in Figure 4. The fea- 
ture does not change in either position or shape with 
thickness. We can exclude first-order or higher-order 
phonons as its origin. Moreover, it is visible only at 
near-resonance excitation (~1.96 eV) which is very close 
to the direct band gap at K points in the Brillouin zone. 
This implies that the feature might be related to vertical 
electronic transition crossing direct band gap. 

First the exciton mechanism could be excluded due to 
the energy scale difference. The so-called Al and Bl 
excitons in bulk M0S2 are characterized with the bind- 
ing energies of ~42 and 134 meV— , and recent optical 
experiments indicated that the two excitonic states ex- 
hibit little change in few-layer M0S2 compared to bulk 
materials^. The exciton binding energies are much 
larger than that of the broad feature around 5 meV. Im- 
purity or defect states are also unlikely to contribute to 
this feature from the energy scale. Impurity and defects 
states appearing in optical measurements are of the same 
energy scale as the exciton binding energj£. It should be 
noted that our samples and those used in the optical 
measurements are from the same source (SPI Supplies). 
The similar conditions of impurities and defects could 
be expected in our samples if existing. Collective elec- 
tric carrier oscillation, so-called plasmon mode and other 
collective modes such as charge density wave (CDW) are 
also excluded from the fact of thickness-independence, 
as Coulomb screening and interlayer coupling would be 
strongly dependent on thickness. 

The potential mechanism responsible for the broad fea- 
ture may lie in the spin-orbit coupling in conduction band 
at K points in their Brillouin zone. Recent band structure 
calculation a 26 ' indicated that spin-orbit coupling con- 
tributes to the splitting at the top of valence band at K 
points [v\ and i/2), which are the origins of Al and Bl ex- 
citons. At the same time, spin-orbit coupling also induces 
a subtle splitting of ~4 meV at the conduction band at 
K points. The subtle splitting would be reflected in a 
resonant-like Raman procedure crossing the direct band 
gap at K points. This naturally explains why it appears 
only at near resonant excitations. The thickness indepen- 
dent Raman feature manifests the thickness-independent 
spin-orbit coupling induced splitting in conduction band. 

In conclusion, two low-frequency modes, the shear 
mode and the compression mode are investigated in the 
present Raman measurements. They exhibit completely 
opposite evolution with the thickness of ultrathin MoS2- 
The pair of low-frequency provides a new and more pre- 
cise method to characterize the thickness of atomically 
thin M0S2 crystal. We also observe a broad feature at 
~5 meV which is likely related to the electronic struc- 
ture. We point out that the feature may originate from 
the spin-orbit coupling induced splitting at the conduc- 
tion band at K points in their Brillouin zone. 
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